Experimental joint signal-idler quasi-distributions and photon-number statistics for 

mesoscopic twin beams 
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Joint signal-idler photoelectron distributions of twin beams containing several tens of photons 
per mode have been measured recently. Exploiting a microscopic quantum theory for joint quasi- 
distributions in parametric down-conversion developed earlier we characterize properties of twin 
beams in terms of quasi-distributions using experimental data. Negative values as well as oscillating 
behaviour in quantum region are characteristic for the subsequently determined joint signal-idler 
quasi-distributions of integrated intensities. Also the conditional and difference photon-number 
distributions are shown to be sub-Poissonian and sub-shot-noise, respectively. 

PACS numbers: 42.50.Ar,42.50.Dv,42.65.Lm 



I. INTRODUCTION 

The process of spontaneous parametric down- 
conversion 0, [3, [I[ is one of the fundamental nonlinear 
quantum processes that can be understood in terms of 
created and annihilated photon pairs. This highly non- 
classical origin of the generated optical fields is responsi- 
ble for their unusual properties. They have occurred to 
be extraordinarily useful in verification of fundamental 
laws of quantum mechanics using tests of Bell inequali- 
ties [3J, generation of Greenberger-Horne-Zeilinger states 
[H , demonstration of quantum teleportation || , quantum 
cryptography [QL dense coding, and many other 'quan- 
tum protocols' Q. Fields composed of photon pairs have 
already found applications, e.g. in quantum cryptogra- 
phy [|| or metrology @. Description of this process has 
been elaborated from several points of view for beams 
containin g ju st one photon pair with a low probability 
S E3) ITll Il2l. Il3l as well as for beams in which many pho- 
ton pairs occur jl4j ]. Also stimulated emission of photon 
pairs has been investigate d flBl [Lsl . \vf\. 

Recent experiments 0, [l9|, HO, HH, HH (and references 
therein) have provided experimental joint signal-idler 
photoelectron distributions of twin beams containing up 
to several thousands of photon pairs. Extremely sensi- 
tive photodiodes, special single-photon avalanche pho- 
todiodes [23| , super-conducting bolometers J23k time- 
multiplexed fiber-optics detection loops [25], l26l. l27l. |28| , 
intensified CCD cameras fl9l . [2^|, or methods measur- 
ing attenuated beams [30l . l3l| are available at present as 
detection devices able to resolve photon numbers. Also 
homodyne detection has been applied to determine in- 
tensity correlations of twin beams [3^, These ad- 
vances in experimental techniques have stimulated the 
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development of a detailed microscopic theory able to de- 
scribe such beams and give an insight into their physi- 
cal properties. A theory based on multi-mode descrip- 
tion of the generated fields has been elaborated both 
for spontaneous (34[ as well as stimulated processes [35| . 
This theory allows one to determine the joint signal-idler 
quasi-distribution of integrated intensities from measured 
joint signal-idler photoelectron distributions. Consider- 
ing phases of multi-mode fields generated in this sponta- 
neous process to be completely random, the joint signal- 
idler quasi-distribution of integrated intensities gives us 
a complete description of the generated twin beams. As 
a consequence of pairwise emission the quasi-distribution 
of integrated intensities has a typical shape and attains 
negative values in some regions. This quasi-distribution 
has been already experimentally reached [HI, H3] for twin 
beams containing up to several tens of photon pairs but 
with mean numbers of photons per mode being just a 
fraction of one. Here, we report on experimental determi- 
nation of the joint signal-idler quasi-distribution of inte- 
grated intensities for twin beams containing several tens 
of photons per mode. Such system may be considered as 
mesoscopic and this makes its properties extraordinarily 
interesting for an investigation. 

Photon-number distributions and quasi-distributions 
of integrated intensities provided by theory are contained 
in Sec. II. Section III is devoted to the analysis of exper- 
imental data. Conclusions are drawn in Sec. IV. 



II. PHOTON-NUMBER DISTRIBUTIONS AND 
QUASI-DISTRIBUTIONS OF INTEGRATED 
INTENSITIES 

In the experiment, whose layout is sketched in Fig. [1] 
(for details, see [HI), the third harmonic field (wave- 
length 349 nm and time duration 4.45 ps) of an ampli- 
fied mode-locked Nd:YLF laser with repetition rate of 
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FIG. 1: (Color online) Sketch of the experimental setup. 
Nd:YLF, amplified ps-pulsed laser source; BBO I, nonlinear 
crystal; /, /i,2, lenses; PH12, 100 fim diameter pin-holes; 
Di,2, pin detectors; BPi^, band-pass filters; ND, adjustable 
neutral-density filter. The boxes on the right-hand side in- 
dicate the parts of the signal amplification and acquisition 
chains. 



500 Hz (High Q Laser Production, Hohenems, Austria) 
is used to pump parametric down-conversion in a BBO 
crystal (Fujian Castech Crystals, Fuzhou, China) cut for 
type-I interaction. The down-converted beams at wave- 
lengths of 632.8 and 778.2 nm are selected by two 100 fim 
diameter apertures and directed into two amplified pin 
photodiodes (S5973-02 and S3883, Hamamatsu Photon- 
ics K.K., Japan) using lenses of appropriate focal lengths 
(see Fig. [1]). The output current pulses are digitized and 
processed by a computer. The overall detection quan- 
tum efficiencies, rj, are 55 % for both arms. First ((to)) 
and second ((to 2 )) moments of photoelectron distribu- 
tions for both signal and idler beams as well as correlation 
of photoelectron numbers in the signal and idler beams 
are obtained experimentally. Additive noise present dur- 
ing detection can be measured separately and subtracted 
from the measured data. The measured moments of pho- 
toelectron numbers can be corrected also for the overall 
quantum efficiency and we then obtain the moments for 
photons. Symbol (m) (("2)) means mean photon number 
in signal (idler) field, {n\) ((n 2 .)) denotes the second mo- 
ment of signal- (idler-) field photon- number distribution, 
and (nin2) gives correlations in the number of signal and 
idler photons. We note that moments of photon- number 
distributions are obtained using the relations: 



(m) = (m t )/r], 

(nf) = (m?)/rf - (1 - r])(mi)/r] 2 , 
(nin 2 ) = (m 1 m 2 )/rj 2 . 



7= 1,2, 



(1) 



Moments of integrated intensities can be directly de- 
rived from moments of photon numbers: 



(Wi) = (m), 

(Wf> = (n?)-(n;>, 7 = 1,2, 
(WiW 2 ) = (nm 2 ). 



(2) 



Multi-mode theory of down-conversion developed in [34j 
using a generalized superposition of signal and noise pro- 
vides the following relations between the above men- 



tioned experimental quantities and quantum noise coef- 
ficients B\, B 2 , D12, and the number M of modes: 



(Wi) 
((AWi) 2 ) 
(AW1AW2) 



MBf, i 
MID12I 2 . 



1,2, 



(3) 



The coefficient Bi gives mean number of photons in mode 
i and D12 characterizes mutual correlations between the 
signal and idler fields. Inverting relations in Eqs. ([3]) we 
arrive at the expressions for parameters B\, B 2 , M, and 
Dia: 



Bi = ((AWi) 2 )/ (Wi), 
Mi = (W 4 ) 2 /((AVF l ) 2 ), 7 = 1,2, 
|Di 2 | = ^(AWiAW 2 )/M. 



(4) 



As follows from Eqs. ([!]), the number M of modes can 
be determined from experimental data measured either 
in the signal or idler field. This means that the experi- 
mental data give two numbers Mi and M 2 of modes as 
a consequence of non-perfect alignment of the setup and 
non-perfect exclusion of noise from the data. On the 
other hand, there occurs only one number M of modes 
(number of degrees of freedom) in the theory (34| as it 
assumes that all pairs of mutually entangled signal- and 
idler-field modes are detected at both detectors. Precise 
fulfillment of this requirement can hardly be reached un- 
der real experimental conditions. However, experimental 
data with Mi M 2 can be obtained. 

Joint signal-idler photon-number distribution p(ni, n 2 ) 
for multi-thermal field with M degrees of freedom and 
composed of photon pairs can be derived in the form 



p{ni,n 2 ) = 



1 (-K) n2 {Bi +K) ni - n2 
r(M) (l + Bi +B 2 + K)^+ M 

g T(m +M + r) 



r— max(0,ri2 — n\ ) 

\{Bi+K)(B 2 +K)Y 
(-K) r (l + Bi + B 2 +K) 



r\(n 2 — r)\(ni — n 2 + r)\ 



(5) 



Determinant K introduced in Eq. ([5]) and given by the 
expression 



K = BiB 2 - \Di 2 \ 2 



(6) 



is crucial for the judgement of classicality of a field. Neg- 
ative values of the determinant K mean that a given 
field cannot be described classically as in case of the con- 
sidered field composed of photon pairs. In Eq. {HJ), the 
quantities Bi+K and B 2 +K cannot be negative and can 
be considered as characteristics of fictitious noise present 
in the signal and idler fields, respectively. The theory for 
an ideal lossless case gives K = —B\ = —B 2 together 
with the joint photon-number distribution p(ni,n 2 ) in 
the form of diagonal Mandel-Rice distribution. On the 
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other hand inclusion of losses and external noise results 
in non-diagonal photon-number distribution p(nx,ri2) as 
a consequence of the fact that not all detected photons 
are paired. We note that pairing of photons leads to 
higher values of elements p(ni, n 2 ) of the joint signal-idler 
photon-number distribution around the diagonal n\ = n 2 
that violate a classical inequality derived in . 

A compound Mandel-Rice formula gives the joint 
signal- idler photon- number distribution p(nx,n 2 ) at the 
border between the classical and nonclassical characters 
of the field (K = 0): 



p{ni,n 2 ) 



r(m + n 2 + M)B^B. 



»>2 

2 



T{M)m\m\ (1 + Bi + B 2 ) ni+n 2 +M ' 



(7) 



If the number M of modes is large compared to mean val- 
ues (ni) and (n 2 ) (i.e. for B\, B 2l and I-D12I being small) 
the expression in Eq. ||7J) can roughly be approximated 
by a product of two Poissonian distributions. If K > or 
K < 0, weak classical or quantum fluctuations remain in 
this Poisson limit of a large number M of modes, as fol- 
lows from the normal generating function [34[ in the form 
G(Ai,A 2 ) «exp(-Ai(ni)-A 2 (n2) + AiA 2 Af| J Di2| 2 ) valid 
in this approximation. Thus there are always mode cor- 
relations in agreement with the third formula in Eq. ([3]) . 

Declination from an ideal diagonal distribution 
p{n\,n 2 ) caused by losses can be characterized us- 
ing conditional idler-field photon-number distribution 
Pc,2(n 2 ; ni) measured under the condition of detected n\ 
signal photons and determined along the formula: 

00 

Pc, 2(n 2 ;ni) = p(nx,n 2 )/ y]p(ni,fc). (8) 

Substitution of Eq. {5| in Eq. ([5]) leads to the conditional 
idler-field photon- number distribution p Ci2 with Fano fac- 
tor F c . 2 : 

F c , 2 {nx) = 1 

, (1 + M/m)[(B 2 + K)/(l + B,)} 2 - {K/Brf 



(1 



• M/m)(B 2 + 
« l + K/B x . 



K)/{l + B 1 )-K/B 1 



(9) 



As an approximate expression for the Fano factor F Ci2 
in Eq. @ (valid for K tzs —B 2 ) indicates nega- 
tive values of the determinant K are necessary to ob- 
serve sub-Poissonian conditional photon-number distri- 
butions. Sub-Poissonian conditional distribution p c 2 
emerges from the formula in Eq. ([5]) that is a sum of 
positive terms in this case. For the ideal lossless case, 
K = —B\ = —B 2 holds and the Fano factor F c2 equals 
0. On the other hand, positive values of the determinant 
K mean that the sum in Eq. ([3J contains large terms 
with alternating sings (this may lead to numerical errors 
in summation) and so the conditional distribution p Cj2 is 
super-Poissonian. For instance, for K small compared to 
Bi, F Cj2 = 1 + (B 2 + K)/(l + Bi). We note that, in this 
approximation, the value of Fano factor F c 2 equals the 



value of coefficient R quantifying sub-shot-noise correla- 
tions and being defined in Eq. (fTI]) below. 

Pairing of photons in the detected signal and idler fields 
leads to narrowing of distribution p_ of the difference 
n\ — n 2 of signal- and idler-field photon numbers: 
00 

P-( n ) = ^2 •Wi-nsK™!)^) (10) 
ni ,7i2=0 

and S denotes Kronecker symbol. If variance of the differ- 
ence ni — n 2 of signal- and idler-field photon numbers is 
less than the sum of mean photon numbers in the signal 
and idler fields we speak about sub-shot-noise correla- 
tions and characterize them by coefficient R [2l| : 



([A(ni - n 2 )] 2 ) 

H — ; r ; r < 1. 

(fll) + (n 2 ) 



(11) 



Joint signal-idler photon-number distribution p(ni, n 2 ) 
and joint signal-idler quasi-distribution Pi(W\, W 2 ) of in- 
tegrated intensities belonging to normally-ordered oper- 
ators are connected through Mandel's detection equation 
111: 



p(ni,n 2 ) 



ni\n 2 \ J 



dW x \ dW 2 W? 1 W 2 



x exp(-Wi - W 2 )Pi(W 1 ,W 2 ). (12) 



Relation in Eq. (|12p can be gen eralized to an arbitrary 
ordering of field operators [3J, |34| | and can be inverted in 
terms of series of Laguerre polynomials. Range of con- 
vergence of these series is determined under the condition 
s < Sth where Sth is given in Eq. (jT5]) later. These series 
define quasi-distributions for s > Sth- 

Provided that an s-ordered determinant K s , K s — 
B ls B 2s - \D 12 \ 2 {B hS - Bi + (1 - s)/2, i = 1,2), is 
positive the s-ordered joint signal-idler quasi-distribution 
Ps(Wi, W 2 ) of integrated intensities exists as an ordinary 
function [34j | which cannot take on negative values: 



P.{Wi,W 3 



1 



T(M)KM V \Di2 1 2 

{B 2s W l /B ls + W 2 )B U 



x exp 



Symbol Im denotes modified Bessel function and T 
stands for T-function. 

If the s-ordered determinant K s is negative, the joint 
signal-idler quasi-distribution P s of integrated intensities 
exists in general as a generalized function that can be 
negative or even have singularities. It can be approxi- 
mated by the following formula (34|: 



P S {W U W 2 ) 



A{W{W 2 ) 



(M-l)/2 



TrT(M)(B ls B 2s 
A 



)M/: 



exp 



2B ls 



W 2 
2B, 



x sine 



^Wj - w 2 



(14) 
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sinc(a;) = sin(x)/x. Oscillating behaviour is typical for 
the quasi-distribution P s written in Eq. (|14p . 

There exists a threshold value Sth of the ordering pa- 
rameter s for given values of parameters Bi, B2, and D±2 
determined by the condition K s = 0: 

s th = 1 + B x + B 2 - ^{B l + B 2 y - IE- (15) 

— 1 < sth < 1- Quasi-distributions P s for s < s t h are or- 
dinary functions with non-negative values whereas those 
for s > s t h are generalized functions with negative values 
and oscillations. 

Similarly as for photon numbers we can define quasi- 
distribution P Si _ of the difference W\ — W2 of signal- and 
idler-field integrated intensities as a quantity useful for 
description of photon pairing: 

/>oo />oo 

P,,-(W) = dWy \ dW 2 
Jo Jo 

x 5(W-Wi +W 2 )P S (W 1 ,W 2 ). (16) 

Quasi-distribution P s _ oscillates and takes on negative 
values as a consequence of pairwise character of the de- 
tected fields if s > s t h- 

There exists relation between variances of the differ- 
ence fii — n 2 of signal- and idler-field photon numbers 
and difference W\ — W 2 of signal- and idler-field inte- 
grated intensities: 

<[A(m - n 2 )f) = ( ni ) + (n 2 ) + ([A(W 1 - W 2 )] 2 ). (17) 

According to Eq. (TIT)) negative values of the quasi- 
distribution P Si _ (as well as these of quasi-distribution 
P s ) are necessary to observe sub-shot- noise correlations 
in signal- and idler-field photon numbers as described by 
the condition R < 1 in Eq. (TTT]) . 

III. EXPERIMENTAL DISTRIBUTIONS 

As an example, we analyze the following experimen- 
tal data appropriate for photons and derived from the 
experimental data for photoelectrons given in [2l| using 
relations in Eqs. {TJ [77 = 0.55]: 

(rai) = 959.21, (n 2 ) = 1078.3, 
(nf) = 971829.7, {nj) = 1218608, 
(mn 2 ) = 1088083. (18) 

These data thus characterize photon fields, as they have 
been obtained after correction for the nonunit detection 
efficiency. Formulas in Eqs. ([2]) and ([4]) then give mean 
number B\ of signal photons per mode, mean number 
B 2 of idler photons per mode, number M\ of signal-field 
modes, and number M 2 of idler-field modes: 

B x = 52.95, B 2 = 50.81, 

Mi = 18.11, M 2 = 21.22. (19) 



Numbers Mi and M 2 of modes given in Eqs. ([T9|) and de- 
termined from data characterizing signal (Mi) and idler 
(M2) fields slightly differ owing to experimental imper- 
fections. That is why we use a mean number M of modes 
[M = (Mi + M2)/2] and determine the coefficient D12 
along the relation in Eqs. Q: 

M = 19.66, \D 12 \ = 52.29. (20) 

Determinant K given in Eq. © then equals -44.23, i.e. 
the measured field is nonclassical. Coefficient R de- 
fined in Eq. (fTTj) equals 0.19 (-7.2 dB reduction of vac- 
uum fluctuations) and this means that fluctuations in 
the difference n\ — n 2 of signal and idler photon num- 
bers are below shot-noise level. This also means [see 
Eq. ([T7J>] that variance ([A(W / i — W^)] 2 ) of the difference 
of signal- and idler-field integrated intensities is nega- 
tive (([A(Wi - W 2 )] 2 ) = -1654). Negative value of this 
variance is caused by pairwise character of the detected 
fields, which leads to strong correlations in integrated 
intensities Wi and W 2 . Also the valu e of covariance C 
(C = (AniAn 2 )/v/([Ani] 2 )([An 2 ] 2 )) of signal n x and 
idler n 2 photon numbers close to one (C = 0.997) is evi- 
dence of a strong pairwise character of the detected fields. 
We note that also a two-mode principal squeeze variance 
A characterizing phase squeezing and related to one pair 
of modes can be determined along the formula: 

A = 1 + Bi +B 2 -2|£» 12 |. (21) 

We arrive at A = 0.18 using our data in Eq. (p?Tj) and so 
the generated field is also phase squeezed. 

The joint signal-idler photon-number distribution 
p(ni,n 2 ) determined along the formula in Eq. ([5]) for 
values of parameters in Eqs. (fT^|) and is shown in 
Fig. O Strong correlations in signal-field ni and idler- 
field n 2 photon numbers are clearly visible. Nonzero ele- 
ments of the joint photon-number distribution p(ni, n 2 ) 
are localized around a line given by the condition m ~ n 2 
as documented in contour plot in Fig. [5] 

Conditional distributions p Cj2 of idler-field photon 
numbers n 2 conditioned by detection of a given num- 
ber ni of signal photons defined in Eq. (JSJ) are also sub- 
Poissonian (see Fig. [3]). The greater the value of the num- 
ber ni of signal photons the smaller the value of Fano fac- 
tor _F Cj 2 given in Eq. If mean numbers {n\} and (n 2 ) 
of signal- and idler-field photons are small compared to 
the number M of modes the joint photon-number distri- 
bution p(ni , n 2 ) behaves like a product of two Poissonian 
distributions and so F c , 2 w 1. Fano factor F Ct2 reaches 
its asymptotic value after certain value of the number ri\ 
of signal- field photons [see discussion below Eq. (0]. 

Strong correlations in signal-field n% and idler-field n 2 
photon numbers lead to sub-Poissonian distribution p_ 
of the difference ri\ — n 2 of photon numbers defined in 
Eq. pi) (see Fig. ED. 

Joint signal- idler quasi-distributions P s (Wi, W 2 ) of in- 
tegrated intensities differ qualitatively according to the 
value of ordering parameter s (sth = 0.15 for the ex- 
perimental data). Nonclassical character of the detected 
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FIG. 2: Joint signal-idler photon-number distribution 
p(ni,n2). 

1.2 n 
1.0 - 

0.8 - 



o.o 1 1 1 1 1 1 . 1 . 1 

100 200 300 400 500 

FIG. 3: Fano factor F c ,2 of the conditional signal-idler photon- 
number distribution p Cj 2 as a function of the number n\ of 
detected signal photons. 



fields is smoothed out (K s = 2.66 > 0) for the value of s 
equal to 0.1 as shown in Fig. EJa). On the other hand, 
the value of s equal to 0.2 is sufficient to observe quantum 
features (K s — —2.53 < 0) in the joint signal-idler quasi- 
distribution P S (W 1 ,W 2 ) that is plotted in Fig. Mb). In 
this case oscillations and negative values occur in the 
graph of the joint quasi-distribution P S (W\, W 2 ). 

Determination of the number M of modes for the over- 
all field has to be done carefully because it might happen 
that the theory shows nonphysical results. There are 
three conditions determining the region with nonclassi- 
cal behavior: K < 0, K + B x > 0, and K + B 2 > 0. 
These conditions can be transformed into the following 
inequalities: 

B^B 2 < \D 12 \ 2 < B 1 B 2 +mm(B 1 ,B 2 ), (22) 

where min denotes minimum value of its arguments. If 
sub-shot-noise reduction in fluctuations of the difference 
m — n 2 of signal- and idler-field photon numbers is as- 
sumed (implying ([A(VFi — VF2)] 2 ) < 0), even stronger 




FIG. 4: Distributions p-(n) of the difference n of signal- 
field (ni) and idler-field (712) photon numbers; n — ri\ — n-2- 
Solid curve without symbols characterizes the experimental 
data. Solid curve with o gives the distribution obtained 
from the joint signal-idler photon-number distribution in the 
form of product of two independent Poissonian distributions 
with mean photon-numbers given by experimental data and 
is shown for comparison. 

conditions can be derived: 

BiB 2 < {B 2 + B 2 )/2< \D 12 \ 2 < BiB 2 + min(Bi,5 2 ) 

(23) 

and we therefore need to fulfill the inequality: (Si — 
B 2 ) 2 < 2 mm(Bi, B 2 ). Assuming B\ > B 2 we arrive at 
the final condition: 

B X <B 2 + yf2B~ 2 . (24) 

The condition in Eq. (|24|) gives limitation to the low- 
est possible physical value of the number M of modes. 
Increasing the value of number M of modes from this 
boundary value the field behaves nonclassically first and 
then its properties become classical. 

Nonclassical character of the detected field is given by 
the condition K < in theory. In experiment we usu- 
ally measure coefficient R given in Eq. (TIT]) in order to 
prove nonclassical character of the field given by the con- 
dition R < 1 . According to the developed theory [34j , if 
the field is classical (K > 0), then there is no sub-shot- 
noise reduction in fluctuations of the difference of signal- 
and idler-field photon numbers (R > 1). On the other 
hand, situation is more complicated for nonclassical fields 
with K < 0. Provided that B\ — B 2 — B, negative 
value of the determinant K implies B 2 — \Di 2 \ 2 < and 
([A(Wi - W 2 )] 2 ) = 2M(B 2 - \D l2 \ 2 ) < 0. Thus the use 
of relation in Eq. (TP?)) gives R < 1, i.e. we have sub- 
shot-noise reduction of fluctuations in the difference of 
photon numbers. If B\ =/= B 2 , it may happen that R > 1, 
i.e. non-classicality of the field is not observed in sub- 
shot-noise reduction of fluctuations of the difference of 
photon numbers. We note that even conditional photon- 
number distributions p c ^ 2 can remain sub-Poissonian in 
this case. 

The above discussion has been devoted to statistical 
properties of photons. Qualitatively similar results can 
be obtained also for photoelectrons. Quantum Burgess 
theorem assures that sub-Poissonian photoelectron dis- 
tribution occurs provided that photon-number distribu- 
tion is also sub-Poissonian [F m — 1 = rj(F n — 1), F n (F m ) 
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a) 



b) 



7-, 
6 

a? , 





means Fano factor for photons (photoelectrons) , r\ is de- 
tection efficiency]. Photoelectron distributions are nois- 
ier compared to photon-number distributions and that 
is why nonclassical properties of photoelectron distribu- 
tions are weaker. 



IV. CONCLUSIONS 



Nonclassical character of mesoscopic twin beams con- 
taining several tens of photon pairs per mode has been 
demonstrated using experimental data. Joint signal- 
idler photon-number distribution, its conditional photon- 
number distributions, distribution of the difference of 
signal- and idler-field photon numbers, and joint signal- 
idler quasi-distributions of integrated intensities have 
been determined to provide evidence of non-classicality 
of the detected twin beams. 
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